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We discuss the coupled variations of the gravitational, strong and electroweak coupling 

O ! 

CNJ , constants and the current knowledge of the nuclear equation of state based on heavy ion 



collision experiments and neutron star mass-radius relationship. In particular we focus in 
our description on phcnomenological parameters, R, relating variations in the quantum chro- 
modynamics scale Aqcd and the fine structure constant a, and 5, relating variation of v, 
C^H, the Higgs vacuum expectation value and the Yukawa couplings, h, in the quark sector. This 

parametrization is valid for any model where gauge coupling unification occurs at some (un- 

^M ' 

specified) high energy scale. From a physically motivated set of equations of state for dense 

matter we obtain the constrained parameter phase space (i?, S) in high density nuclear envi- 
ronments. This procedure is complementary to (although currently less powerful than) those 

> ■ .... . 

^\ • used in low density conditions. For variations of Aa/a = 0.005 we find that the obtained 

~f\ ' constrained parameter lies on a strip region in the (i?, S) plane that partially overlaps some 



of the allowed values of parameters derived from primordial abundances. This may be of 
interest in the context of unification scenarios where a dense phase of the universe may have 



existed at early times. 
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I. INTRODUCTION 

The properties of matter are the result of the interplay of fundamental interactions. Current 
cnowledge points towards new physics beyond the standard paradigm of three fundamental forces 
l|], and suggests that unification should occur at high — but so far unknown — energies. In that 
context, possible variations in the strong, electroweak and gravitational interactions are not inde- 
pendent of each other but, in a wide class of unification scenarios they are coupled. In this sense 
a variety of limits on the variation of the fine structure constant, a, and the ratio of quark masses 
to the quantum chromodynamics (QCD) scale, m q / 'Aqcd, have been established. In particular 
for the latter one there are recent measurements from consideration of Big Bang nucleosynthesis 
(BBN), quasar absorption spectra, and the Sm isotopes in the Oklo natural nuclear reactor which 
was active about 1.8 billion years ago (for a review see e.g. 2j). 

In order to model possibly coupled variations of fundamental constants in a framework valid 
for any model where gauge coupling unification occurs at some (unspecified) high energy scale 
it is assumed that the electroweak scale is derived by dimensional transmutation and that the 
varying couplings are due to some dilaton-type scalar field. In this way a set of phenomenological 
parameters (R, S) |3J relating Yukawa couplings, h, Higgs vacuum expectation value, v, QCD 
mass scale, Aqcd, and fine structure constant, a is usually assumed. More in detail, S is defined 
as Av/v = SAh/h and since there is some model dependence a safe assumption is allowing S 
to variate. Typically S ~ At where At is the sensitivity from the top quark and defined as 
At = d In mw/<9 In a,\ , where mw is the mass of the W boson and <Zj are the input parameters of the 
supersynnnetnc model following the Bartueri-Giudice measure for an observable flEJ. Assuming a 
value for the Higgs vacuum expectation value of v = Mp exp(— 8-7r 2 c/h^ ) [6J], where c and ht central 
values can be considered to be of order unity c ~ ho ~ 1 and Mp is the Planck mass, and variating 
this expression it results that Av/v ~ I6OA/1//1 which will be taken as a reference for the canonical 
choice S ~ 160 from now on. Considering supersymmetric models with more massive particles may 
enlarge S up to ~ 500. As for the R parameter in some grand unification models it is predicted 
that Aqcd, is modified as 5(m/ Aqcd)/ (m/ 'Aqcd) ~ Rbci/a where m is the quark or electron 



mass in t 



v 



lese models. This result is strongly model dependent and it could yield R ~ 36 for some 



However, the large coefficients in these 



models [7[ or even negative sign values as argued in 
expressions are generic for grand unification models and therefore signal that possible variation 
of masses and the strong interaction scale may be easier to detect than the variation of the fine 
structure constant a. 



Previous works [3( [9( focused on the low baryonic mass density environment of BBN where 
Pb ~ 10~ 5 g/cm 3 and physical quantities of interest such as the nucleon mass difference in vacuum 
Q = ro „ - ,n p , the neutron lifetime and the h,nd,„ g energy of deuteriun, For e ra mp>e, Coc et 
al. [3|] found that assuming a dilaton model and coupled variations one could accommodate the 
observational data of primordial abundances of light nuclei with an scenario of coupled variations 
for several parameter choices, including Ah/h = 1.5 10 -5 , R = 16, S = 240 and Ah/h = 2.5 10~ 5 , 
R = 45, S = 240; in both cases Aa/a = 2Ah/h. 

In this work we use the complementarity provided by the extreme conditions of matter described 
by the nuclear equation of state (EoS) in heavy ion (HI) and neutron star (NS) physics. Even if 
the precision of inferred nuclear observables differs from the BBN nuclear abundance data or other 
techniques [l(| , one can assess the constraining power of these high density environments and scan 
the (a, R, S) parameter space. A subsequent goal (which is beyond the scope of this paper) will 
be to carefully study particular implications for unification scenarios. 

In an astrophysical scenario the EoS of nuclear matter is crucial for describing NS structure 
and, in particular, to obtain its mass-radius relationship. Conversely, kinematical observables, 
photospheric radius expansion or thermal emission measurements for NSs allow to infer compact 
masses and radii that constrain the possible values for central densities and the EoS itself. Typical 
NSs have on average a mass M^s ~ 1-4Mq and radius .Rats' ~ 12 km. Roughly, their structure 
can be described as an external crust where matter has densities below saturation density, no ~ 
0.145 fin - or, correspondingly, in mass density po ~ 210 14 g/cm 3 , and an internal core where 
densities are larger and so far unknown. These objects are born in the aftermath of a supernova 
event with birth energies AE ~ E grav ~ 5 „ NS ~ 10 53 erg released when a progenitor exhausts 
the core reactions fueling the star and it can no longer support its own gravity. Some of them may 
have extreme external magnetic fields (of the order B ~ 10 9 — 10 15 G) and emit regular radiation 
pulses with periods T ~ 10~ 3 — 10 1 s, being known as pulsars. The interior structure of NSs in the 
spherically symmetric static approximation is obtained by solving the Tolman (TOV) equations 



111 ], given an EoS P = P(e,T) relating pressure, P, energy density, e, and temperature T. 

This same EoS is also relevant for understanding the properties of nuclei and yields of HI 
collisions (at higher temperatures) experimentally accessible on Earth. Both sides of low and high 
densities allow to partially test the phase space of matter in the density-temperature plane. 

In this work we consider the physics of matter under the extreme conditions that may resemble 
somewhat those of the early universe in the nuclear dense phase, by studying how the EoS can be 
used to constrain simultaneous coupled variations of various dimensionless fundamental couplings. 



Such a possibility, in particular a time variation of the fine structure constant, steams from ob- 
servations of quasar absorption systems and a fit of positions of absorption lines |12| suggesting a 
smaller value in the past. Additionally, there are some indications that current data from Keck 
and VLT could be reinterpreted as a possible spatial variation or gradient in the fine structure 
constant a [la . Il4j |. In this sense some of the parameters in models already describing coupled 
variations are obtained in the context of primordial abundances in BBN low conditions, these probe 
different conditions that for the high density interiors of NSs and HI collisions. As we explain in 
this work the possible variation of masses of different population species under beta equilibrium, 
namely nucleons and electrons, can be modified by the non vanishing values of a variation in a. 
The in-medium effects can be parametrized by a meson-field model of strong interaction, due to the 
fact that high density environment is now considered, and typically one has to consider "dressed" 
or effective masses. In this way gravitational, weak and strong interactions are considered into the 
picture at high density conditions at the phenomenological level. The modelization used has to 
deal with a number of uncertainties but the conclusions remain generic otherwise. 

In section [TT] we introduce the phenomenological model used in this work to describe the set 
of EoSs for nuclear data and the generalized formalism of in-medium variations of quantities of 
interest in the high density conditions. In section IIIII we analyze the constraints of the (R, S) 
parameter phase space to matter conditions as described by the set of EoSs used in this work. We 
apply our constraining procedure and discuss the agreement or overlap of some of the obtained 
values of (R, S) parameters to those previously quoted in the literature. Finally, in section IIV} we 
summarize and give some conclusions. 

II. VARYING COUPLINGS 



Previous works trying to assess the importance of varying fundamental constants in the stellar 
structure [15l . 1 161 ] have considered a polytropic EoS P = Kp , where p is the mass density, V = 
1 + 1/n and n is the polytropic index. Matter at high densities can be usually considered as a 
degenerate system where T ~ 0, since Fermi energies Epi of the degenerate it/i-particle species 
are much higher than thermal energies Ept/hBT 3> 1. The interplay of the constituent particle 
interactions is in-built in the EoS and further influences other observables as in HI collisions [171 ] or 
in a compact object mass and radius configuration [18(. Typical central densities for these systems 



are not accurately known but for purely hadronic stars it is assumed 1914211] that they could reach 



up to « 5no- About the composition, it has been hypothesized that more exotic objects composed 



of strongly interacting nucleonic, hyperonic or deconfined quark matter at their central regions [22] 
may also exist. 

There are recent indications, coming from spectroscopic measurements along the line of sight 
of quasars, of space and time variations of a 13J]. These time variations are at the parts-per- 



million level. However, one must realize they apply to very low density environments. The only 
plausibly realistic way to explain such spatial variations would be with a chameleon-type field, 
and in that case the variations will also be environment-dependent, -they will depend on the local 
density 23, [2J]. Therefore, the values of the couplings at high densities attained in the center of 
NSs or in HI collisions can differ from those on Earth or on low-density environments (as probed 
by quasar absorption systems or the cosmic microwave background (CMB)). In other words, the 
parts-per-million results of |13J ] need not apply to the high density contrasts in NSs, and in what 
follows we will consider variations Aa/a ~ ±10 -3 . Those can be accomodated by the CMB data 

n 

in an analysis performed by Menegoni et al. [25J who obtained a result —0.013 < Aa/a < 0.015 at 
95% C.L. from WMAP 5-year date combined with ACBAR, QUAD and BICEP experiments data 
or in the framework of some other models [26[][lJ]. Again we stress that, as pointed out in [14i. |23|]. 
there is no inconsistency with a value of a differing in such a proportion as those compared to 
quasar absorption since they could probe different conditions. 

We now address the variation of fundamental constants in terms of density dependent couplings. 
To parametrize varying couplings we will use the fact that when a dimensionless coupling such as 
the fine-structure constant is changed by a small amount a = ao(l+5 a ), changes in other quantities 
can be related to it through coefficients AX /X = kxAa/a for X = Xo(l + kx8 a ). This is the case 
in unification scenarios shown in [3J, for which these changes can be phenomenologically described 
by two parameters: R, relating AAqqd and Act, and S relating Av and Ah (the Yukawa couplings, 
all assumed to be the same) that we will discuss below. 

Matter in a NS is subject to gravitational, electroweak and strong interaction. Therefore, we 
will consider coupled variations in the fine structure constant a = e 2 /he and, additionally, allow 
for variation in particle masses on the lepton and hadron sectors. We define a set of dimensionless 
quantities to work with, namely the proton to electron mass ratio [i = m p /m e , and for physically 
meaningful changes in the gravitational constant G we consider the ratio ac = Gmj/hc = rrif/Mp 
for the ith-type particle species. From the work of [3] the variations in \x can be written as 
Afi/fi = [0.8R - 0.3(1 + 5)] Aa/a. 

On physical grounds both R and S may be expected to be positive, however there are some 



simplified phenomenological descriptions 



26[ where one has, for example, S = — 1 and R = 109 



6 



(although they have some drawbacks as compared with quasar data). Let us emphasize that there 
is no experimental data which unambiguously shows what the correct unification model is like (or 
indeed, if one wants to be skeptical, whether high-energy unification happens at all). In this work 
we have adopted a phenomenological approach, treating R and S as free parameters, and ask at 
what extent they are constrained by current knowledge of nuclear EoS physics, and how robust are 
these constraints. A list of models, representative of several possible unification scenarios, where R 
can have positive or negative values of several hundreds, are discussed in 



27, 



281 ] . As mentioned 

before, in the case of S its value is related to the neutralino mass and may be as large as S ~ 500 
as pointed in 29(. Here we conservatively assume that R, S can take values, either positive or 
negative between the maximum absolute value, R,S£ [—500, 500]. As most observables of interest 
will depend on a linear combination of R and S we phenomenologically assume the same range of 
variation for both, as this will allow us to more easily identify possible degeneracy directions. 
Despite some earlier results, there is currently no similar evidence for varying /i from a detailed 



analysis [30 



311 ] . If these results are correct they indicate that A/i//i <C Aa/a and therefore (in 
our class of models) R ~ 3(1 + S)/8, which could be used to eliminate one parameter from the 
analysis; however we will proceed without this constraint. 



Let us notice that there is a vast number of model EoSs in the literature 



321 ]. In this work 



we will firstly describe highly degenerate charge-neutral isospin asymmetric relativistic matter 



by a non- linear Walecka model (NLWM) 33j including baryons (neutrons and protons), mesons 



(Lorentz scalar, vector and iso- vector, a, oj and p respectively) and leptons (electrons). We choose 



a parametrization, TM1 



341 ] . which considers isospin symmetry since we have verified that more 



advanced parameterizations as PK1 



351 ] yield same trend results. Since a full comprehensive 



treatment is not affordable due to the large amount or current treatments of these systems we have 
accounted for the spread of the nuclear EoS by using a phenomenological approach. In section IIIII 
we define a pressure spread parameter, 5 and allow it to vary from S = to 5 = 0.5 to partially 
size the constraining power of this kind of description based on effective field theories. Let us point 
out that the main goal of each EoS is to describe the nuclear and astrophysical phenomenology 
currently known and there is no consensus about the actual degrees of freedom that should be 
present. Conservatively, we consider the NLWM a suitable EoS since it has proven to give a 
reasonable description of nuclear phenomenology but we keep in mind that the interior of NSs is 
largely unknown. In [3|] they use a non-relativistic potential model where the binding energy of 
deuterium is based on the nucleon and a and oj meson mass. We introduce additionally the p meson 
in our treatment since our system is largely neutron rich and therefore non isospin symmetric. 



The variation in the meson-nucleon field coupling in the NLWM, <?j, (i = a,oj,p) can be 
parametrized as gj = g^mf/Mp), so that Agi/gi = Ami/mi. The relative variations in 



the electron and proton masses are given by 



3(, Am e /m e = 0.5(1 + S)/Aa/a, Am p /m p 



[0.8-R + 0.2(1 + S)] Aa/a, and since the nucleon mass difference parameter is Q = m n — m p 
and variations are AQ/Q = (0.1 + 0.75 — 0.6R)Aa/a we find that for neutrons Am n /m n = 
[(0.1 + 0.75 - 0.6R) + (m. p /m n )(0.1 - 0.55 + 1.4R)] Aa/a. Note that in the isospin-symmetric 
case we recover Am n /m n = Am p /m p ; moreover, requiring that the mass difference still vanishes 
as the couplings vary leads to the consistency condition 6-R = 75 + 1. 

In the dense medium, the interacting Fermi systems of baryons have effective masses m* = 
m p ~ daP, m n = m n ~ 9v a anc ^ the relative mass variation can be obtained for protons as, 



(0.8^ + 0.2(1 + 5)) — -^ 
a g a 

and in an analogous way for neutrons. 



Am* p 
m* 



m* g a 



To include the meson mass variation we consider the Feynman-Hellmann theorem [36| where the 
expected value of the quark condensate (in the light sector u, d, s) for the proton state is obtained 
as (p\qq\p) = dm p /dm q . Using a description of the a meson as a SU(3) singlet [9( the meson mass 
variation can be related to the s-quark mass, m s , as dm a /dm s = (a\qq\a) = 2/3. Taking values 
for the quark masses (using h = c = 1) m u ss 4 MeV, m^ w 6 MeV, (m u « m,j ss m q ), m s sa 104 
MeV we obtain, accordingly, for the a meson: 

Am a _ 2 m s Am s 2 {m u + m d ) Am, 



for the u) meson, 



m a 3 m a m s 3 m, 



Am,,, m s A?n s 



and for the p meson 



rricj m^j m s 



Am p m s Am s 



(3) 



(4) 



m p m p m s 

Using the assumption that all relative variations of the Yukawa couplings are similar [3j we 
have Am Stq /m Stq = 0.5(1 + S)Aa/a. Assuming small variations of the dimensionless couplings 
in this dense system, we now describe the set of equations to solve for a charge neutral spin- 
saturated beta-equilibrated system given an input baryonic particle number density, n. The self- 
consistent relativistic mean fields values a = (a), uj° = (uj°), p° = (p°) are obtained from the set 



a = g a {l + Ag a /g a )n s /m'l , lo° = ^(1 + Ag^/g^n/m'l, m 2 p p° = g p (n p - n n ) /2, fj, n = /j, p + p e . 
Particle number densities of different involved species are rij = {n p ,n n ,n e } and each one is given 
in terms of the Fermi momentum, kpi as rij = k Fi /3ir 2 . Conserved baryonic number densities and 
electrical charge neutrality imply additionally n = n p + n n , n p = n^. The effective meson masses, 



m 'o, m 'uj, include non- linear self-interaction terms in the NLWM [33]. 

The explicit expressions for the chemical potentials are given by p p = E F + g^oo® + g p p° /2, 
H n = Ep + g^u — g p p°/2, p e = E F where E F , E F and Ep are, respectively, the proton, neutron 



and electron Fermi energies given by E F = \/k Fi + m* 2 (l + 2 m J ), and m* = m e . Finally, the 
scalar density is given by, 

ii*. < 1 ^) t ' a (5) 

If the solution exists, it is given by values solving the self-consistent set of non- linear equations 



for the NLWM EoS P = P(e, T = 0). In this way the pressure is given by [33J] 



1 ^ f k ™ k A dk 1 v^ /■*« A; 4 rffc 



«=p,n 



where the mesonic field pressure is given in terms of the non-linear couplings k, A and £ as, 



-1 o ,. .Ara„, o 1 _ .Am,,n, n^ 1 o,. .Aro„ w n ,o 1 i 1. ^ 1 



p ™ = ^< 1+2 ^f >" 2 +?"" (1 + 2 ^7><^ +?" 2 < 1 + 2 ^r ) W-a»"-^+^(^. 



and the energy density is, 



* = ^ e r* ^ w* 2 + m?2(i + 2 ^ } + ^ ^ r e dk ^ 2 + m?a + 2 ^ } + em ' 

i=p,n V l «=e 

(8) 
where the mesonic contribution is e TO = —P m - 



III. RESULTS 

We now discuss the results obtained from exploring the (a, R, S) space of the dense system. We 
note that our analysis differs from that of Coc et al. [3] where some unification models were studied, 
by choosing fixed values of a, R and S. Here our goal is to scan the phenomenological (R, S) space 
for fixed choices of a, in order to ascertain the feasibility of constraining this space using the current 
knowledge of the nuclear EoS. A fine-structure constant variation of Aa/a = +0.005 is considered. 



In Fig {T] we show the variation of particle population for protons (dashed line) AY P or neutrons AY n 
(solid line) defined with respect to the unchanged value of a, Al$ = Y. a ' a — Y i a ' a ~ i = p, n as a 
function of particle number density. We use Aa/a = +0.005 and a 'canonical' choice of parameters 
R = 20,5 = 160 [3J]. It can be seen that the particle population changes in an appreciable way 
decreasing the number of protons as a is enhanced. In a similar fashion to [14| a change in a or 
quark masses translates into a change in BBN abundances. 
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FIG. 1. Variation of particle population fractions with Aa/a 
baryonic number density. 



-0.005, R=20, S=160 as a function of 



In order to illustrate the current experimental constraints on the nuclear EoS we plot in Fig. [2] 
the pressure versus the baryonic particle number density from a variety of input data. With solid 
line we plot EoS data from the work of Steiner et al. [18[ deduced from a selected group of NS mass 
and radius measurements. EoS data from HI collisions has also been plotted from the symmetric 
nuclear matter (SNM) case (dashed line) and pure neutron matter (PNM) case (dotted line). Note 



that NS matter is in beta equilibrium and therefore none of the two idea 
system close to the PNM case. Both data seem to suggest a soft EoS |37|, 



cases apply, but it is a 



38]. 



There is a vast number of parameterizations of the nuclear EoS in the literature, including several 
treatments of the nuclear interaction. From those using the non-relativistic potential interaction 
models 



39l | to the Fermi gas or the relativistic fields as Muller and Serot (£ = 0, £ = 0) [40j or the 



Shen j34l | TM1 EoS. They all aim to describe a many-body system to give insight to the nuclear 



observables and the astrophysically deduced masses and radii of NSs. In Fig. [3] (upper panel) we 
include a representation of some of the most representative nuclear EoS as obtained for the SNM 
and PNM cases and for beta equilibrium TM1 for the sake of comparison. We also plot the EoS 
NS constraints (dotted line). We see that most EoS do not describe the whole range of densities 
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FIG. 2. Experimental nuclear EoS constraints in pressure from NS mass-radii measurements 
collisions in the SNM and PNM case [17| as a function of baryonic number density. 



and HI 



in complete agreement. Some of them are not appropriate as e. g. the PNM EoS of Muller and 
Serot (C = 0,£ = 0), since it is too stiff or the Fermi gas EoS since it does not provide enough 
pressure at large densities. We can see in the lower panel in Fig. [3] again a representation of the 
same set of EoS but using the HI constraints (dashed boxes) for the SNM and PNM cases. We see 
that again the TM1 EoS allows the representation the neutron rich side of those boxes, however 
some of the Muller and Fermi gas EoS do not describe the tendency of these regions as happened 
with NS constraints. 

In order to partially consider the spread of the nuclear EoS we introduce a phenomenological 
parameter 5 to parameterize the softness of the EoS and its effect is shown in Fig. [H In this way for 
the range < 6 < 0.5 the pressure is modified accordingly as P(l — 6). Therefore 5 = corresponds 
to no modification of the pressure as given by the EoS considered and 5 = 0.5 a 50% softened EoS. 
Note that since our reference EoS describes better the stiff side of the constraints, we have 5 > 0. 
We consider the phenomenological spread in the nuclear EoS and see that the unchanged TM1 EoS 
describes the neutron rich systems and as 5 > more isospin symmetric systems can be accessed. 
With this procedure we partially size the combined effect of the uncertainty in the isospin content 
and the softness of the EoS and the robustness of our findings. 

Having gained a feeling for the effect of variations in a, we now explore the R, S parameter 
space. Interestingly, and although the effects are relatively mild for the 'canonical' choice of R and 
S, we will now see that there are regions of parameter space where there is no possible solution to 
the set of equations showing, therefore, 'exclusion regions'. 

In Fig. [5] we show for a value Aa/a = 0.005 and R, S in the interval [—500, 500] using TM1 
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FIG. 3. Some frequently used nuclear EoS and NS constraints as deduced from 
constraints as deduced from 17j (lower panel). 
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FIG. 4. Set of phenomenologically softened TM1 EoS (S = 0,0.2,0.5) versus EoS coming from HI and NS 
constraints. 
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the maximum and minimum values for pressure as a function of density for selected values ti/uq = 
1, 2, 3, 4, 5 and for spreads 5 = 0, 0.2, 0.5, crosses, circles and triangles respectively. The aim of this 
figure is to show the pressure spread and that not all unification scenarios, i.e. values of R and S 
are consistent simultaneously with existing EoS deduced from NS and HI datasets in the bounded 
areas. This is due to the sensitivity of the nuclear EoS to the strong effect on the effective particle 
masses in the dense medium from varying the couplings. 
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FIG. 5. Maximum and minimum values of pressure spread in the phcnomcnologically softened TM1 (5 = 
0, 0.2, 0.5) when R, S are varied in the [—500, 500] interval for Aa/a = 0.005 as compared to EoS from NS 
and HI constraints. 



In Fig. [6] we plot with crosses, circles and triangles values of R and S yielding pressure values 
compatible with EoS from NS constraints using TM1 EoS with Aa/a = 0.005 and applying a 
spread of 5 = 0, 0.2, 0.5 respectively. We see that, globally, the values compatible with EoS from 
NS constraints lie on a variable width strip. Without being fully comprehensive on the scan of 
the allowed R and S values we can see that equally stiff or soft EoS have constraining power in 
their phase space. Some of these values follow the same robust tendency when an extreme 5 = 0.5 
variation in the EoS is considered and they are superimposed. However one must note that the 
experimental uncertainty of these data is much larger that in the low density BBN case. Note 
for example that there are quantities that are especially sensitive to the parameter change as, for 
example, the Li abundance that is sensitive to a relative change of 50% when a is changed 1% 
[141 ] . The same constraints are applied in the HI collision data show in Fig. [7] We see that in 
either case the tendency is the same, showing robustness. The range of R shown in both figures is 
somewhat lower than the originally prescribed due to exclusion of those larger values since they do 
not provide a valid self-consistent solution of the set of population equations, namely due to the 
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contribution of effective masses of the particles. 
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FIG. 6. R, S phase space exclusion region. Crosses, circles and triangles denote values of pressure compatible 
with NS constraints as results from applying a spread of 5 = 0, 0.2, 0.5 for Aa/a = 0.5%. 



Note that as the density increases the uncertainty of these measurements is larger and there is 
a trivial constraint on R, S values since pressures must remain positive. We note that NS matter is 
mostly neutron matter but indeed not pure neutron nor symmetric matter (those cases constrained 
in [17[ or at low density in e.g. [43J]) but we use these constraints as an indication of the behavior. 

Finally following the analysis performed in [3J we find that both set of parameters Ah/h = 
1.5 10" 5 , R = 16, S = 240 and Ah/h = 2.5 10" 5 , R = 45, S = 240 with Aa/a = 2Ah/h are 
indeed contained in our set of possible solutions. Using their assumed value Aa/a ~ 10~ 5 would 
result in pressure deviations smaller than 1%, and therefore it would produce a negligible variation 
on the EoS. Therefore the constraining power of the procedure presented in this work is meaningful 
for values of Aa/a larger than this and currently present in some theoretical models of unification 



and data analysis 



25] 



As mentioned, we expect that some mild dependence on the EoS modelization may arise, in 
particular a more enriched description could somewhat enlarge or reduce those constrained regions 
but would not produce a change of tendency in the constrained values since these descriptions 
use the same basic strong interaction, quark masses, fine structure and gravitational constant 
ingredients. 
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FIG. 7. Same as in Fig. ©but for HI PNM (top panel) SNM (bottom panel) case for Aa/a = 0.5%. See 
details in the text. 

IV. CONCLUSIONS 



In this work we have performed a study of combined variations of the gravitational, strong 
and electroweak coupling constants based on the current constraints on the nuclear high density 
EoS. Using a phenomenological approach and assuming a variation of the fine structure constant 
compatible, for example, with recent analysis of CMB data, Aa/a ~ 10~ 3 , we find that the valid 
solutions provided a variation of couplings in the studied model, lie on a variable width strip in 
the phase space of R, relating variations in the Aqqd an d the fine structure constant a, and S, 
relating variation of v, the Higgs vacuum expectation value and the Yukawa couplings, h, in the 
quark sector. If we further restrict our results to positive R, S values on the first quadrant, based 
on more stringent physical assumptions for the most feasible grand unification models, we find 
that some of the parameters already obtained in the context of primordial abundances in BBN 
low conditions remain valid when considered in the high density conditions of NS interiors and 
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HI collisions. Although the experimental precision of nuclear observables in HI collision and NS 
physics is lower than in nuclear finite systems or primordial nuclear abundances, the extended 
system EoS depends crucially on in-medium density effects. Therefore we find that EoS from HI 
and NS data may constrain the variation of a, R, S parameters at densities beyond uq and provide 
additional non-trivial constraints on them. 

Our results show the potential for using the high density region of the density phase space of 
matter to partially constrain fundamental physics in addition to low density tests using the CMB. 
Future work to explore the relevant parameter space in more detail, and discuss how these tests 
can complement other constraints from stellar reactions is needed. 
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